Rate theory for correlated processes: Double-jumps in adatom diffusion. 
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We study the rate of activated motion over multiple bar- 
riers, in particular the correlated double-jump of an adatom 
diffusing on a missing- row reconstructed Platinum (110) sur- 
face. We develop a Transition Path Theory, showing that the 
activation energy is given by the minimum-energy trajectory 
which succeeds in the double-jump. We explicitly calculate 
this trajectory within an effective-medium molecular dynam- 
ics simulation. A cusp in the acceptance region leads to a yT 
prefactor for the activated rate of double-jumps. Theory and 
numerical results agree. 



Reaction rates and diffusion rates in crystalline en- 
vironments typically have an Arrhenius dependence on 
temperature r\ — z/exp (—Ets/T). This asymptotic rate 
at low temperatures may be calculated using Transition 
State Theory (TST) § , where E TS is the energy of 
the saddle-point atomic configuration separating the ini- 
tial and final states (the Transition State), and v is a 
temperature-independent prefactor involving the curva- 
tures of the energy surface. 

Many rates and transitions are not described by a sim- 
ple crossing of a single barrier, and we should expect 
that their rates will not be given by the simple Arrhe- 
nius form. In this paper we study a double-jump: a cor- 
related diffusion event where an atom crosses two barri- 
ers. The development of field-ion microscopy (FIM) and 
scanning tunneling microscopy (STM) has made it possi- 
ble to track the motion of individual atoms at surfaces Q , 
and to directly measure the rates of these correlated tran- 
sitions [j3|,|]]. Developing what we call Transition Path 
Theory (TPT), we show the rate is determined by the 
energy Etp of the transition path (minimizing the en- 
ergy among all paths which succeed in the double-jump), 
and has the asymptotic form T2 — C\/Texp (—Etp/T). 
We develop an efficient numerical method to calculate 
this minimum-energy path and the rate, and use it to de- 
scribe double-jumps along the troughs of the missing-row 
reconstructed Pt(llO) surface within effective-medium 
theory [pi — making contact with the Arrhenius double- 
jump rate measured in the recent Arhus experiment Q. 

Previously, the rates of double-jumps have been dis- 
cussed theoretically as so-called dynamical corrections to 
transition state theory J6j,(jj or in the context of general- 
ized Langevin equation models 



FIG. 1. The geometry of the Pt(110)-surface with an 
adatom and the Transition Path (TP) for a double-jump. 
At upper left, one sees the actual volume of atoms used in 
the simulations, with periodic boundary conditions and static 
atoms at the bottom. Lower left shows the geometry of the 
system with the adatom #1 in a valley formed by the recon- 
struction. At right we show schematically the TP and the 
energy. The TP starts at t = — oo with the system locally 
(close to the adatom) at the transition state (TS) potential 
energy saddle-point between wells A and B; in addition, there 
is an energy AE = Etp — Ets stored in the degrees of free- 
dom far away. The additional energy radiates in from infinity 
and helps the adatom cross the well and get to the second TS. 
Our calculation shows that the adatom rolls off the hill, hits 
atoms 2&4 and knocks them aside, then is hit from behind 
by atoms 3&5, who boost the adatom up the other hill. As 
t — > oo the energy AE is again radiated away to infinity jToj . 

Consider an extra adatom on the reconstructed (110) 
surface (Fig. |IJ). A single-jump diffusion event is a ther- 
mal fluctuation in which the adatom moves along the 
troughs in the surface from one potential well over the 
transition barrier to a neighboring well (say from well 
A to B in the right-hand Fig. [I]) . In a double-jump, 
the atom moves from one potential well over two barriers 
before settling down in a well (A to C in Fig. 0). The 
double jumps can be distinguished from two subsequent 
single-jumps only if there is a separation of time scales: 
the duration time for a double-jump to pass through the 
central well (roughly a picosecond for Pt) has to be much 
smaller than the inverse rate 1/ri for single jumps (1.4 
nanoseconds at the highest temperatures we work at). 

In the case of a single-jump the transition state en- 
ergy Ets is the smallest possible thermal fluctuation for 
a single-jump to occur. For our calculation, E^s is 469.5 
meV. What is the smallest possible thermal fluctuation, 
in which the adatom performs a double-jump: fluctuat- 
ing to the top of one hill, sliding down, up to the top 
of a second hill, and then down? Here the energy of the 
minimal fluctuation involves both kinetic and potential 
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energy: we must look in phase space for the transition 
path (TP): the lowest energy trajectory which will give 
us a double jump. 

Before we more formally explain how the TP is cal- 
culated we shall describe the qualitative features of the 
path. Consider first the time evolution of the system if 
it is started out in a TS and the adatom is given a slight 
push. The adatom will then slowly leave the saddlepoint 
and slide down the hill into a well but it will not make it 
over the next barrier because energy will be transferred 
to other atoms in the system. Eventually the adatom will 
be at rest at the bottom of the well and all the energy 
£/rs will be in the degrees of freedom far away from the 
adatom. The TP is defined so that it brings the adatom 
to the top of the next barrier (the next TS) with the 
minimum additional energy (AE — E^p — -^ts) neces- 
sary. In the TP the system is at t = — oo at the TS for 
all the degrees of freedom close to the adatom, but an 
energy AE is stored in the degrees of freedom far away. 
As the adatom slides down the hill the energy is radiat- 
ing in from infinity helping the adatom up to the second 
TS (Fig. |l|). As t — > oo the adatom (and all local co- 
ordinates/momenta) approaches the second TS and the 
energy AE is again radiated away to infinity pj. The 
additional energy AE results in a higher activation en- 
ergy for double-jumps than for single-jumps. 

In finding the minimum energy double-jump path, it is 
important to use as variables the positions and momenta 
r at t = when the adatom is crossing the bottom of the 
central well, because at this time the additional energy 
AE is most localized near the adatom. Let p x be the 
x-momentum of the adatom at t — 0, where x is the 
coordinate along the trough. To find Ttp of the TP, 
we specify the positions and the momenta of the nearby 
atoms within a radius R of the adatom, and place atoms 
outside R at their relaxed positions. We then vary p x 
to find p™ m , where the adatom just barely succeeds in 
the double-jump: crosses into the well on the right at 
large times, and into the well on the left at large negative 
times [0. We calculate the total energy, and minimize 
with respect to the positions and momenta of the nearby 
atoms. With all atoms in their relaxed positions, 652.0 
meV of energy is needed in p x to get a double-jump. 
Optimizing degrees of freedom within R = 1, 2.2, 3.1, 
3.8, and 4.9 nearest neighbor distances of the adatom, the 
total energy decreases to 621.7, 595.8, 591.4, 590.2 and 
589.8 meV respectively. The energy converges rapidly as 
more degrees of freedom are optimized and T — > Ttp- We 
estimate Etp to be 589.8 meV: 120.3 meV higher than 
the single-jump Ets- The calculated TP appears to have 
both of the obvious possible symmetries: one reflection 
plane, and one reflection plane plus time-reversal. The 
time-reversal symmetry implies that the local atoms both 
start and end at zero velocity at the saddle points. 

To calculate the rate at low temperatures, we need 
to sum over all the low-energy double-jump trajectories. 
Fig. ^ gives an idea about how the double-jump region 
in phase space looks. The figure shows this region in 



a 2-dimensional cut in phase space (the plane spanned 
by the x-coordinate of a bridge atom (atom 2 in Fig. |l|) 
and p x of the adatom; all the remaining coordinates are 
at their Ttp values). R is the region corresponding to 
trajectories where the adatom will jump over the barrier 
to the right; L is the region where the adatom came over 
the barrier from the left. The TP is the lowest-energy 
trajectory which came from the left and makes it to the 
right, and is situated on the cusp at the (co-dimension 
two) intersection of the bounding surfaces of R and L. 
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FIG. 2. Energy contours (100 meV) in a cut in phase space, 
nip,, and dm are the atomic mass and nearest neighbor dis- 
tance of Platinum. The bridge atom is atom 2 in Fig. 1. All 
other coordinates are fixed at their values in the Transition 
Path (TP) at time t — 0. L(R) indicates regions in phase 
space, where the adatom makes a jump to the left(right), and 
the solid lines are trajectories where the local atoms end at the 
saddlepoints. The double-jump region (gray shaded) is above 
both lines, and so the bounding surface for the double-jump 
region has a cusp passing through Ttp , the phase-space point 
taken by the TP at t = 0. 

Our method for calculating the rate of double-jumps 
is analogous to the standard TST method for calculating 
the rates for single-jumps. TST expresses the hopping 
rate as the flux through the dividing surface; at low tem- 
peratures it is well described as a harmonic expansion of 
the energy about the TS 
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In this expression (. . .) denotes a thermal average, Zyy 
denotes the partition function in one well. The reaction 
coordinate x r equals x\ at the top of the barrier; the w's 
denote the eigenfrequencies at the bottom of the well, 
and the w TS 's denote the eigenfrequencies found by the 
harmonic expansion at the TS (excluding the imaginary, 
unstable frequency of the saddle point). Choosing the 
coordinate x r to be along the eigendirection of the imag- 
inary frequency minimizes the recrossing corrections jjj . 
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Similarly, the double-jump rate can be calculated as 
the thermal average of the flux through a surface x = 
due to double-jump trajectories: 



r-2 = (Qd 3 v x S(x))/Z w . 



(2) 



Here Qdj is one if the trajectory is a double-jump, and 
zero otherwise. In both cases the rate is calculated as a 
flux through a surface, with a 9-function keeping only 
trajectories which succeed in the transition. 

We generate an importance-sampling distribution of M 
trajectories, {r Q }, in the neighborhood of the transition 
path by adding Ttp to a thermal ensemble restricted to 
x = at the bottom of the well. We found it important 
numerically to choose i to be a relative coordinate: the 
difference between the x-coordinate of the adatom and 
the average x-coordinate of the four bridge atoms (#2-5 
on Fig. |l|). We then calculate the ratio of the double- 
jump to single-jump rates: 
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We determine which trajectories are double-jumps nu- 
merically: we run the molecular dynamics trajectory for- 
ward and backward in time for a few picoseconds until 
the trajectory either recrosses the original well (failure, 
Qdj = 0) or crosses the bottom of both of the adjacent 
wells (success, Qdj = !)• This expression is complicated 
by the fact that we are measuring the flux for single- 
jumps at the top of the barrier, and the flux for double- 
jumps at the bottom of the well: the term in square 
brackets is precisely the ratio of probabilities of being 
at these two planes, in the harmonic approximation (the 
u)'s being the eigenfrequencies in the plane x — 0, and we 
evaluate Yl^jJ S /u>i to be 1.722). Our double-jump term 
is fully nonlinear. 

Fig. U shows the "TP MD" values calculated using 
equation |3|. Also shown are rate calculations performed 
using the more traditional method described by Voter 
and Doll Q (transition state molecular dynamics "TS 
MD"), where trajectories are started with a thermal dis- 
tribution at the TST dividing surface, and the fraction 
of double-jumps is directly measured. Within the sta- 
tistical error bars the two methods give the same rate, 
as they should. However, the uncertainty in the tradi- 
tional rate determination increases drastically at lower 
temperatures because the fraction of double-jump tra- 
jectories in the thermal ensemble becomes small. In our 
language, the old method centers attention not at Ttp, 
but far away at the single-jump transition state. 
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FIG. 3. Arrhenius plot of the rate of double-jumps rel- 
ative to single jumps, showing the Transition State Molecu- 
lar Dynamics (TS MD), Transition Path Molecular Dynamics 
(TP MD), and Transition Path Theory results. The TP The- 
ory curve is r%lr\ = C2VT exp — (Etp — Ets)/1cT, where the 
one parameter C2 is fitted to the TP MD data to be 0.022 
meV" 1/2 . 

We now discuss the asymptotic behavior of the double- 
jump rates at low temperatures. In analogy to the har- 
monic expansion around the TS which lead to the TST 
rate with an Arrhenius behavior for the single-jump rate 
(Eq. [j]) we can define a transition path theory rate r2 PT 
for the double-jumps by performing an expansion of 
equation || around the TP. 

Qdj is a Heaviside step function in p x — p™ m (T), 
(Fig. ||) , where f is the phase space point at t = exclud- 
ing x and p x of the adatom. The x and p x parts of the 
phase-space integral in the numerator of Eq. || (choosing 
here for convenience x = x) can be explicitly carried out, 
leading to 



T2 = ^- I dr ex P (-£(r)/T). 



(4) 



where E(f) = E(x = 0,p™ in (f), f ) is the energy on 
the boundary of the double-jump region. In the low- 
temperature limit a harmonic expansion of E(T) around 
the TP in the 6A^ — 2 variables T can be performed and 
the resulting multidimensional Gaussian integral can in 
principle be carried out. Focusing on the temperature de- 
pendence of the prefactor, we note that p™", and hence 
E has a cusp in one degree of freedom (see Fig. ||), and 
that there is a 6N — 3-dimensional subspace, where the 
energy is presumably locally quadratic. These 6N — 3 de- 
grees of freedom should give rise to a factor VT (each) in 
the prefactor, whereas the cusp degree of freedom gives 
by integration a factor T, so overall we get the prefactor 
T x T^ N -^' 2 x T/T 3N = VT. The TPT rate for the 
double-jumps is therefore of the form 



TPT 
r 2 



dVTeM-ETp/T), 



(5) 



with C\ a temperature-independent constant. Alterna- 
tively we may write rJ PT /r TST = C2VT exp(— (Etp — 
Ets)/T), where C2 is the constant C\jv. 
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The TPT-rate expression is shown as the solid curve in 
Fig. H with only the constant C2 fitted to the calculated 
rates. The TPT rate is clearly in agreement with the 
simulation results even up to the highest temperature 
T = 65 meV. At high temperatures corrections caused 
by anharmonic effects can be expected and higher order 
multiple jumps (triple, quadruple, etc.) may also play a 
role in the diffusivity. It is easy to fit the TP MD data of 
Fig. U with an Arrhenius form with constant prefactor, 
but it yields the wrong energy barrier (139 meV, 16% too 
high). 
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FIG. 4. Temperature dependence of TP MD data for the 
rate of double-jumps relative to single-jumps. The black dots 
show the data corrected for the predicted exponential tem- 
perature dependence (f(T) = AE/T, AE = E T p - E T s 
= 120.3 meV). The curvature of this set of points proves 
that the prefactor depends on temperature, and is well de- 
scribed by the solid line, which is a plot of I11C2VT, where 
C 2 is 0.022 meV~ 1/2 as in Fig. 3. The open triangles 
shows the data now also corrected for the \/T-dependence 
(/(T) = AE/T - lnT/2). TP Theory predicts this set 
of points this to be a constant (I11C2), consistent with 
the data. Linear regression (dashed line) gives a slope of 
+0.47 ± 0.2 meV ~ 0.4% of AE. This slope is a reasonable 
estimate of how much further our estimated Etp could be 
lowered by further optimization. 

Since we sample the vicinity of the TP in our numeri- 
cal method, we can extend the calculation of the double- 
jump rate to much lower temperatures to confirm the 
theoretical asymptotic behavior. (The double-jump rate 
becomes very small.) In Fig. ^ we show the result for 
T%jr\ dividing out first the predicted Arrhenius depen- 
dence (solid symbols), and further the predicted vT- 
dependence (open symbols), ending up with a set of data 
points consistent with a constant. Hence our numerical 
method confirm the TPT rate as given by Equation |[ 

Comparing with the recent experiments |i) on the 
Pt/Pt(110)-system we note that in our model the barrier 
for single jumps is E^s = 469.5 meV while the activa- 
tion energy experimentally is found to be around 0.8 eV. 
It is well known that the EMT potential tends to un- 
derestimate diffusion barriers for Pt [0|. However, the 
TPT analysis confirms the experimental observation of a 
thermally activated form; the calculated additional acti- 



vation energy for double-jumps Etp — Sts = 120 meV 
is in good agreement with the experimentally determined 
value which is of the order 0.1 eV. 

We finally note that we have applied the TPT to the 
traditional Langevin equation (Ohmic damping) previ- 
ously used to describe double-jumps M. Here the TP 
becomes the optimal time evolution of the external "fluc- 
tuating" force. The additional activation energy AE for 
double-jumps equals the work done by friction in the adi- 
abatic potential only in the low friction limit. We observe 
the same y/T prefactor. 
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